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AN INVESTIGATION OF THE FEASIBILITY OF 
OBTAINING LATERAL STABILITY DERIVATIVES 
FOR A LINEAR AIRCRAFT BY MATCHING ANALOG 
COMPUTER TRANSIENT RESPONSES TO FLIGHT 

TEST DATA 

SUMMARY 

Using a North American NAVION as a flight test vehicle, and an 
electronic analog computer as a simulator, a study was made of the 
feasibility of obtaining certain lateral stability derivatives by matching 
the transient responses of an analog computer on which the lateral 
operational equations of motion for a NAVION aircraft were set up, to 
the actual airplane lateral response in sideslip angle, roll rate, and 
yaw rate. 

The matching process was achieved by supplying time histories 
of the actual recorded control motions as forcing function inputs to the 
computer, and systematically varying those coefficients of the analog 
which corresponded to the stability derivatives of the dynamic equations 
of motion. 

The method was demonstrated to be a practical procedure for 
obtaining the lateral stability derivatives of a linear aircraft. The use 
of four different types of control deflections as forcing functions pro- 
vided a means for matching of computer responses which corresponded 
to a wide range of flight maneuvers. 
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The method is considered suitable for further investigation, and 
development for use under different flight conditions and with non- 
linear aircraft. 
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INTRODUCTION 



Since the close of World War II, and the advent of the high speed 
and altitude characteristics of jet aircraft, the need for a quantitative 
knowledge of the dynamics of flight has resulted in rapid development 
of the field of aircraft dynamics. The past six or seven years, in 
particular, have brought considerable progress in the field of dynamic 
stability flight testing of aircraft and missiles. The early theoretical 
and experimental investigations in this area by the Cornell Aeronau- 
tical Laboratory, and by the Instrumentation Laboratory of the Massa- 
chusetts Institute of Technology, have been supplemented on many 
fronts, and in particular, by the NACA, which in recent years has 
developed its dynamic stability testing techniques to the point where 
free-flight testing of rocket powered and telemeter equipped models 
has become commonplace. 

The role of full scale dynamic flight testing is twofold. First, 
it provides a designer with the conventional data on the characteristics 
of an aircraft for the particular conditions of the test. Secondly, and 
in the long view, more importantly, dynamic flight testing provides a 
means for solution of the inverse problem of dynamics, wherein 
the equations describing a system are assumed known, its output in 
response to suitable forcing functions is measured, and thus the co- 
efficients of the equations, or stability derivatives, can be determined. 
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Once established, these stability derivatives are of considerable value. 
They are indispensable to the aerodynamicist who must calculate sta- 
bility characteristics over a wide range of flight conditions, and who is 
thus provided with a measure of the accuracy of the theoretical and 
empirical derivatives used in preliminary design calculations. Deter- 
mination of the stability derivatives from flight tests also aids the 
designer by indicating lines along which the stability may be improved, 
either by geometric design changes, or by artificial augmentation. In 
the latter regard, the successful design of any system for automatic 
control of an aircraft is greatly dependent upon an accurate knowledge 
of the stability derivatives of the airplane. 

The basic approaches to the problem of the determination of 
stability derivatives from dynamic flight testing, and in particular, 
from experimental time histories, can be separated into two categories 
- the "equations of motion" method, which essentially solves the in- 
verse problem described above, and a second category, "response 
curve fitting," which identifies a collection of methods using various 

mathematical and graphical concepts to extract from the data certain 

* 

combinations of stability derivatives and inertia terms called transfer 
function coefficients. All of the individual methods which make up the 
above two categories will not be described in this report, but an excel- 
lent discussion of each current method can be found in Ref. 1, Table I 
presents a summary and comparison of these various methods for 
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analysis of dynamic systems. As can be seen from Table I, the choice 
of a method for analyzing dynamic response data will depend primarily 
on whether transfer functions, or actual stability derivatives are 
desired. 

Several investigations have been reported (ie. Ref. 2) which 
have attempted to obtain stability derivatives by applying a particular 
combination of the above referenced methods which essentially involved 
the conversion of time response to frequency response, followed by a 
matching process which utilized both dynamic response templates and 
an analog computer to produce transfer function coefficients. These 
constants were then shown to uniquely determine the time response of 
an analog of the test airplane which accurately duplicated the time his- 
tory of the actual test vehicle. Evaluation of the actual aerodynamic 
stability derivatives, however, required supplemental calculations, and 
was possible only in the case of those derivatives which exerted the 
strongest influence on the airplane response. 

The major limitations of the above method are the length of time 
required for conversion of the flight data from time response to fre- 
quency response form, and the necessity for the subsequent use of tem- 
plates to determine approximate transfer function coefficients. It was 
considered that both of these limitations could be avoided in the applica- 
tion of the above method if the matching process were to be applied 
directly to the transient flight data. 
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This report, then is concerned with a study of the feasibility of 
obtaining certain lateral stability derivatives by matching the transient 
responses from an analog computer on which the operational equations 
for a NAVION aircraft were set up, to the actual airplane lateral re- 
sponse in sideslip angle, roll rate, and yaw rate. The matching pro- 
cess was achieved by supplying time histories of the actual recorded 
control motions as inputs to the computer, and systematically varying 
those coefficients of the analog which corresponded to the stability 
derivatives of the dynamic equations of motion. 

An additional objective of this investigation was to determine the 
suitability of a pulse width telemeter system as a data collecting 
scheme for a dynamic flight test program requiring very accurate 
determination of the recorded quantities. 

This investigation was carried out at the James Forrestal Research 
Center, Princeton University, Princeton, New Jersey, during the 1956-57 



academic year. 
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EQUIPMENT AND PROCEDURE 



EQUIPMENT 

Airborne Instrumentation 

Test Airplane. The test airplane was a North American NAVION, 
which is an all metal, four place, low wing, single engine airplane. 

The NAVION is equipped with a Continental E-185 engine, rated at 185 
HP at 2300 RPM and 29 inches manifold pressure, at sea level. The 
wing is unswept, and the tail configuration conventional. The control 
surfaces are all of beaded skin construction, and have no trailing edge 
extensions. Aileron and rudder trim tabs are of the fixed bend type. 
The elevator tab is of conventional design, and can be adjusted from the 
cockpit. 

The test airplane was modified by the removal of the two rear 
cockpit seats to provide space for the test instrumentation. (Fig. 6). 
The test airplane had a take-off gross weight of 2870 pounds with full 
service, pilot and co-pilot, and all test equipment. 

A three view drawing and photograph of the NAVION are shown 
as Fig. 1 and Fig. 2, and the general specifications of the test airplane 
are tabulated in Table II. 

Power Supply. Power available for instrumentation in the test 
airplane included both 12 and 28 volt direct current, supplied by a 28 
volt engine driven generator, and two 12 volt batteries connected in 
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series. A 500 volt-ampere inverter, with both frequency and voltage 
regulation, was installed behind the co-pilot's seat to provide 115 volt 
400 cycle three phase alternating current power to the rate gyros and 
the telemeter transmitter. 

Considerable effort was expended in obtaining an inverter which 
could be frequency controlled, since the accuracy of the rate gyros was 
dependent upon the input frequency being the same as the calibration 
frequency of 400 cycles per second. 

A circuit diagram of the aircraft 28 volt system and the inverter 
system is included in Fig. 3. 

Aircraft Instrumentation. Standard aircraft instruments were 
used to determine indicated airspeed and altitude. Measurement of the 
following quantities was required with respect to stability axes: side- 
slip angle, yaw rate, roll rate, and aileron and rudder position angles. 

In choosing the type of components needed to measure the above quan- 
tities, an attempt was made to utilize instrumentation compatible with 
the telemetering equipment. For this reason, the transducers selected 
were of the potentiometer type, not only because they were capable of 
providing signals of a sufficient magnitude, but also because their excita- 
tion source readily provided zero and full scale reference voltages. 

In view of some past history in this type of flight testing of various 
amounts of coupling between the response of the aircraft and certain of 
the instrumentation components, considerable attention was given to the 
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choice of transducers with satisfactory frequency and phase response 
characteristics. These remarks particularly apply to the choice of a 
suitable device for recording sideslip angle. 

In this investigation, sideslip angle was sensed by a light vane 
mounted on a boom approximately one maximum fuselage diameter 
(four feet) ahead of the left wing tip. The moment of inertia of the vane 
was very low, and no gearing was introduced in the linkage in order to 
keep the natural frequency and damping within acceptable limits. An 
Electro-Mechanical Laboratory ultra-low torque potentiometer, of only 
.008 inch-ounces friction torque, was obtained for use with this side- 
slip vane. In Appendix A is presented a detailed analysis of both the 
theoretical and experimentally obtained natural frequency and damping 
ratio of the sideslip vane. 

The sideslip vane had originally been designed to operate through 
an autosyn transmitter, and thus required some modification in order to 
accommodate a potentiometer. 

In order to reduce boom resonance, as had been experience in 
previous tests with the NAVION, the boom for the sideslip vane was 
stiffened by three angle strips equally spaced around the periphery, and 
running the full length of the boom. The sideslip vane and boom are 
shown in the assembled condition in Fig. 4. 

Rudder and aileron angles were measured by 10, 000 ohm Model 
5301 single turn Helipot precision potentiometers. The rudder pot 
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was attached to the rudder itself, and was actuated by a cable under 
spring tension. At the time the instrumentation system was designed, 
the flight test program was planned to include control inputs from the 
rudder only. For this reason, the instrumentation installed for re- 
cording aileron angle was intended only to indicate if the ailerons 
remained locked throughout the rudder input and the resulting tran- 
sient oscillation. Accordingly, the pickoff of aileron angle was made 
from control cable movement, with the pot attached close to the right 
aileron in order to avoid the effects of cable stretch. As is discussed 
in Appendix C, the later decision to include aileron deflections as 
forcing functions required subsequent correction of the recorded aileron 
data. 

The measurement of yaw and roll rates was made by type 36128 
Giannini rate gyros. This type of gyro operates on a power supply of 
115 volts 400 cycle three phase AC to the motor, and 115 volt single 
phase AC to the heater. The transducer is of the potentiometer type, 
with a DC resistance of 5250 ohms. The gyros were installed as near 
the aircraft center of gravity as possible. However, since the gyros 
were of the angular velocity rather than acceleration type, there was 
not the necessity for precise location that is associated with the latter 
type gyro. A level and inclinometer were used to align the gyros so 
that the input axes would be parallel to the appropriate body axis of the 
aircraft. Since deflection of the gyro spring-restrained gimbal was 
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only two degrees at full range, cross -coupling effects were considered 
to be negligible. However, in order to minimize any error due to rate 
coupling, the spin axis of each gyro was aligned parallel to the body 
Y -axis. The gyros were bolted directly to a rigid fuselage member. 
Since the engine vibration frequency was sufficiently above the natural 
frequency of the gyros (17 cycles per second), there was no engine 
noise introduced into the gyro signal. These gyros were obtained on 
loan from the Instrumentation Laboratory of the U. S. Naval Air Test 
Center, NAS Patuxent River, Maryland. 

Signal Circuit Design. The component of the instrumentation 
which controlled the design of the signal circuit was the sideslip vane. 
The requirements for recording the sideslip angle without introducing 
phase differences were negligible torque, and a relatively high natural 
frequency of the vane. The transducer selected has been shown to have 
met these specifications (Appendix A), but the fact that it was a 360 
degree single turn potentiometer pickup, whereas the sideslip vane was 
expected to operate through only a moderate arc, necessitated the use 
of 45 volt DC excitation of the sideslip pot rather than the zero to five 
volts range supplied to the remaining transducers. Since the resolution 
of the particular sideslip pot used was exceptionally good, the accuracy 
of the sideslip system was comparable to that of a system utilizing 
autosyn pickups, but had none of the problems of AC signal circuits 
and was completely compatible with the telemetering and recording 



10 



equipment. The use of a 45 volt DC signal circuit supply, however, 
necessitated a method of supplying this voltage to the sideslip circuit, 
and at the same time, limiting the full scale voltage in the remaining 
signal circuits to five volts. 

In achieving the above, it was desirable to utilize a common DC 
voltage source for all transducer pickups, since the telemeter was to 
sample the particular supply voltage available 20 times per second, 
and correct this to full scale by a continuous calibration. 

La addition to the above requirements on the signal circuit design, 
it was found necessary to bias the sideslip voltage pot in order that the 
sideslip signal voltage would also vary between zero and a full scale 
value of approximately five volts. 

The particular signal circuit design employed is shown in Fig. 5, 
in schematic form. The essential element of the circuit is the bias 
supply voltage potentiometer, pot "A". With the sideslip vane against 
the limit stop in one direction, pot B A" was adjusted so that the output 
voltage of the sideslip pot was zero. The vane was then deflected to the 
opposite limit stop, and the output sideslip pot voltage recorded. Pot 
B B" was then adjusted to bring the full scale reference voltage, and 
thus also the voltage in the remaining signal circuits, in agreement with 
the full scale sideslip voltage. This procedure was repeated two or 
three times, as necessary, in order to nullify any loading on the system. 

Once adjusted, however, pots "A" and "B" were locked, and required 
no further adjustment for the remainder of the project. 
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During the initial flight and transmission tests, considerable AC 
ripple was found to enter the DC signal voltage circuits, very likely 
from the 400 cycle inverter. Condensers C and D in Fig. 5 
were installed, and were found to reduce the ripple to a negligible value 
without otherwise altering the form or amplitude of the signals. 

Fig. 6 shows the general location and arrangement of the air- 
borne instrumentation components. 

Telemeter System. The telemeter unit was of the pulse width 
type, and was designed and built by the Applied Science Corporation of 
Princeton (ASCOP). This unit was the key element in the instrumenta- 
of this project. The basic functional components of this pulse width 
system are transducers, all of the potentiometer type; a time division 
multiplexer, or commutator; and electronic coder, or keyer to put the 
intelligence into pulse form; an RF transmitter; and ground station 
equipment for receiving, decoding, and recording the data in usable 
form. 

The airborne component of the system is designated the ASCOP 
D Series PW Multicoder, and the ground station equipment the ASCOP 
M Series PW Ground Station. A schematic diagram of the functions 
of the pulse width data system, as adapted for use with this particular 
project, is presented in Fig. 7. 

A unique feature and very great advantage of pulse width data sys- 
tems is that due to the multiplexing or time sharing operation, all data 
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channels traverse the same path during transmission. Since zero and 
full scale references are transmitted with the data, this results in a 
calibration signal, effective for all channels, being transmitted 20 
times each second. These signals supply a continuous two-point cali- 
bration, and thus allow for the removal of almost all of the system 
errors, with the exception of errors introduced in the instrument pick- 
ups and the mechanical switch, and they make possible the use of this 
telemeter in tests of long duration and varying environmental conditions. 

The output voltages from the aircraft instrumentation pickups, 
and the calibration voltages described above, are sampled in sequence 
by a rotary switch. At the normal sampling rate of 20 samples per 
second on each of the 45 channels, the data pulse length for any one 
channel is short in comparison to the interval from sample to sample. 
This permits the use of a storage circuit to create a nearly continuous 
voltage output . The output of the sampling switch is a series of data 
pulses whose voltage amplitude is proportional to the amplitude of the 
physical quantity measured by the corresponding pickup. These samples 
are supplied to an electronic coder circuit (keyer) which converts the 
amplitude pulses into pulses having time durations proportional to the 
measured quantities. These pulses are used to produce modulation of 
the transmitted RF carrier. The RF signal is detected at the ground 
station, with the output of the receiver being a train of pulses identical 
in nature to the output of the converter in the airborne equipment. 
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These pulses are separated by a switch running synchronously with 
the switch at the transmitter, with the pulses representing individual 
data channels being placed on separate outputs. The pulses are then 
demodulated to yield continuous voltages representing the intelligence 
transmitted. 

Ground Station Instrumentation 

The ground station equipment included the Telemeter Ground 
Station (described above), a Magnetic Tape Recorder, a GEDA 
Analog Computer, and suitable recorders. These components are 
illustrated in Fig. 8. 

Magnetic Tape Recorder. An Ampex Model 309 C dual track 
three speed tape recorder was an integral component of the ground 
station equipment, and provided a convenient means for recording, 
storage, and playback of the standard pulse width modulated data. 
Provisions were also included for placing voice transmissions directly 
on the magnetic tape from either or both the aircraft and the ground 
station. 

GEDA Analog Computer. The analog computer used in this 
investigation was a Goodyear Aircraft Corporation Model L3 (GEDA) 
linear electronic differential analyzer. Twenty -four automatically 
stabilized DC computing amplifiers were available with open-loop DC 
gain of greater than 5 x 10^ , and negligible drift. The computer 
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incorporated an automatic error indicator, and had a guaranteed 
accuracy of one percent. Provisions were available for accurately 
setting computer board potentiometers by the use of a special null 
indicator , 

Data Recorder. Qualitative analysis of the flight data during 
tape playbacks was made on a Brush two channel recorder. Once the 
particular data run was selected for further analysis, a Sanborn Model 
154-100B four channel recorder was used. This instrument had a very 
high natural frequency (42 cps), and the low frequency response was 
flat to zero cps. 

Instrument Calibration 

Five components of the instrumentation system required calibra- 
tion. These were the two rate gyros, the two control deflection poten- 
tiometers, and the sideslip vane potentiometer. The aircraft altitude 
and airspeed indicating systems were not calibrated, as previous cali- 
brations indicated that the errors in these systems were negligible for 
the conditions of this investigation. 

The two rate gyros were calibrated on a pendulum which consisted 
of a 3/4 inch diameter aluminum shaft and a 12 inch diameter instru- 
ment pan. The pendulum, the overall length of which was 121.75 inches, 
was pinned to the overhead truss structure of the Forrestal hangar. 

A scale was mounted on a small table in such a way that displacement 
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of the pan from its rest position could be measured. The gyros were 
mounted one at a time on a wooden base plate which was clamped to 
the pendulum pan, and their potentiometer outputs were connected to 
a Brush Pen Recorder. The pendulum was displaced from its rest 
position and allowed to swing freely. Because of friction in the system, 
swing amplitudes were measured on two successive swings in the same 
direction and averaged. The two maximum voltages recorded during 
the swing were also averaged and the maximum rate attained during 
the swing was determined from the following relation : 

_ d 31.0 
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where 



P - natural period of the pendulum, determined from the 
Brush Recorder trace, seconds. 

- length of the pendulum, inches. 

d - average maximum displacement of the pan from rest 
position, inches. 



The pendulum was allowed to swing until once again at rest with 
readings taken as described above over the entire range of gyro rates. 
The angular rate gradient of both gyros was found to be linear and equal 
to 0. 610 degrees per second, per percent full scale potentiometer 
voltage. The gyro calibration curves are presented as Fig. 9. It was 
found that the gyros were very sensitive to the frequency of the 400 
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cycle supply voltage and as a result, great care was taken during 
calibration, and also during the actual flight test program, to main- 
tain the frequency of the supply voltage at exactly 400 cycles per 
second. 

The control deflection potentiometers were calibrated by meas- 
uring the resistance between the output terminals over the full range 
of control deflections. The actual control deflections were measured 
with a vernier protractor capable of accuracy to 0.1 degree. The 
rudder deflection gradient was found to be linear and equal to 0.405 
degrees per percent full scale potentiometer voltage. The aileron pot 
calibration was based on total aileron angle, and was performed over 
the entire range of aileron travel in order to determine the extent of 
aileron differential for the NAVION aircraft. The calibration was 
then converted to average aileron angle, and was found to be linear at 
a constant value of 0.1933 degrees per percent full scale, throughout 
the deflection range used in the flight test, which included -10 to 10 
degrees. The amount of aileron differential gearing in this range was 
found to be extremely small. 

The sideslip vane potentiometer was calibrated in the same 
manner as the control deflection potentiometers, and the calibration 
curve is presented in Fig. 10. The sideslip gradient was found to be 
linear and equal to 0.468 degrees per persent full scale potentiometer 
voltage. 



Since the telemeter equipment corrected its output to percent 



full scale voltage and since full scale voltage on all the signal circuit 
potentiometers varied by the same increment as the reference voltage 
for small fluctuations of the reference, it was unnecessary to cali- 
brate the potentiometers through the telemeter equipment. 

Careful consideration was given to the possible errors which 
might be introduced in the data by the dynamics of the recording in- 
struments. The natural frequency of the rate gyros was 17 cycles 
per second and their damping ratio was between .45 and .7. The 
natural frequency and damping ratio for the sideslip vane were experi- 
mentally determined in flight to be 15.44 cps and .081 respectively. 
(See Appendix A.) Since the natural frequency of the aircraft was 
approximately 0.4 cycles per second, the attenuation of the signals 
due to the dynamics of the instruments was effectively zero. 



PROCEDURE 



General 

Aircraft transient responses in the lateral modes of motion were 
obtained by introducing rudder and aileron control motions of the pulse, 
doublet, and step-function types. Rudder pulses, rudder doublets, and 
aileron doublets were found to be the most effective in exciting the 
* Dutch Roll" motion, which in turn produced the most satisfactory 
transient time histories for the investigation in question. Aileron step- 
functions were also used, however, primarily to provide the time 
history of roll-rate response to lateral control motion. 

Ideally, sufficient transient response data for all types of forcing 
functions could be obtained in one flight of approximately an hour’s 
duration. However, due primarily to the need of the investigators to 
gain some proficiency with the instrumentation and technique of applying 
satisfactory forcing functions, several flights were made in the process 
of collecting data. It is worthy of note that the basic scheme of instru- 
mentation and data collection proved singularly reliable during the flight 
test - data collection phase. After initial installation the instrumentation 
required only a very minimum of modification and adjustment, and in 
each instance the trouble was minor in nature and easily remedied. 
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Aircraft 

Preflight preparation consisted of checking the full scale voltage 
applied to the sideslip potentiometer to insure that it matched that of 
the other signal potentiometers, and to insure that the full scale signal 
voltage did not exceed the 4-6 volt range which the telemetering 
equipment needed for automatic scale calibration. The fuel tank was 
checked and if necessary topped-off to insure that the aircraft weight 
at the beginning of each flight was constant. When deemed advisable, 
a complete ground check of all communication and signal transmission 
circuits was made prior to take-off. After completing the ground 
checks all test equipment and the aircraft generator were secured to 
insure maximum engine power for take-off. 

During the climb to test altitude, the generator and inverter were 
energized and power was applied to the rate gyro heaters and to the air- 
borne telemeter heater and filament circuits. Upon reaching test alti- 
tude, voice contact was established with the telemeter ground station 
and all instrumentation circuits were energized. Because of an appar- 
ent deficiency in the inverter frequency control, the inverter frequency 
was monitored continuously and maintained at 400 cps by varying the 
DC voltage supplied by the aircraft generator to the inverter. Voltage 
control was obtained by means of an instrument -panel mounted potenti- 
ometer which was tied into the aircraft voltage regulator circuit (Oper- 
ation of the synchronous motors in the rate gyros at exactly 400 cps was 
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necessary to insure that the output signals followed the calibration 
curves obtained at that frequency). 

All flight tests were conducted at a density altitude of 5000 feet 
and at a true airspeed of 129.2 mph (120 mph IAS at 5000 density 
altitude). A large plot of free air temperature versus pressure alti- 
tude was carried in the cockpit to facilitate adjusting pressure altitude 
(altimeter set at 29.92 in. Hg) with temperature to obtain and main- 
tain 5000 feet density altitude. 

The general technique employed in the actual test runs fits the 
following description. A condition of steady-state flight was established 
at test altitude and airspeed and held for at least 30 seconds prior to 
each control input. This period of steady-state flight insured that the 
aircraft was as nearly trimmed and established in steady state condi- 
tions as possible, and provided a sufficient period of recorded steady 
state signals to allow for circuit zero-signal biasing when the data was 
played back into the analog computer or pen recorder. At the end of 
the steady state period a radio transmission identifying the type of con- 
trol displacement to be made was sent and recorded on the magnetic 
tape. Simultaneously, all controls were locked mechanically in the 
cockpit with the exception of the control to be displaced. Five seconds 
after the above transmission the control forcing function was applied 
and the control was then returned to the original steady-state position 
and locked. All controls were held locked until the resulting aircraft 
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oscillation had essentially subsided. In the case of the aileron step 
functions the displacement was applied as rapidly as possible and the 
control was then held as rigidly as possible. The resulting motion 
in this case was permitted to proceed for 4 to 6 seconds before 
recovery was initiated. 

To lock the elevator and ailerons when forcing the aircraft with 
rudder displacements, a chain and turnbuckle system attached to the 
deck on the right side of the cockpit was used. The copilot's control 
wheel was removed and replaced by a sprocket wheel with a set screw 
to lock it to the control shaft. When the aircraft was established in 
steady-state flight, the sprocket wheel was pushed forward on the con- 
trol shaft until it butted against the instrument panel. The chain was 
then put on the sprocket and the turnbuckle s taken up, making the chain 
taut. Just prior to initiating the rudder deflection, the set screw was 
tightened, locking the sprocket fitting to the control shaft. With the 
chain on and tight, the ailerons were effectively locked, and, with the 
sprocket fitting against the instrument panel, no forward motion of 
the control shaft was possible. The copilot then held forward pressure 
on the sprocket fitting to prevent aft movement of the control shaft. 

To lock the rudder when putting in aileron displacements, and to 
return the rudder to its steady-state position after rudder displace- 
ments, a chain and turnbuckle system on the pilot's right rudder pedal 
was employed. This was attached to the right hand seat track on the 
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pilot te side. Because the aircraft required a slight amount of left 
rudder for trim at the test flight condition, this attachment was also 
useful in establishing trimmed flight by adjusting the turnbuckle. 

Ground Station 

In order to visually monitor the data runs, an external patch- 
board provided on the telemeter receiver (See Fig. 8) was wired in 
such a way as to utilize a voltmeter -type translator for each of the 
recorded variables, and also to allow the decoded output to be placed 
on a pen recorder. The telemeter receiver and the magnetic tape 
recorder had been integrated prior to this investigation, thus requiring 
only the proper positioning of controls on the telemeter console to 
accomplish tape recording of the flight data. Each quantity measured 
during the flight tests was transmitted on three separate telemeter 
channels. Thus, each measured variable was sampled 60 times each 
second. Once the airborne equipment was operating, the telemeter 
receiver was tuned for maximum signal reception, and the "zero" and 
"full-scale" reference voltages were adjusted for proper calibration of 
output signals. 

All pertinent comm uni cations between the ground station and the 
aircraft were recorded on magnetic tape along with the flight test data. 
Prior to commencing the recording of test data on each flight, the 
ground station operator pre-recorded on the magnetic tape the data 
necessary to identify the particular flight, and test conditions as 
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regarded aircraft weight, lift coefficient, and the telemeter channels 
being used for each measured quantity. 

Concurrent with the reception and tape recording of the flight 
test data, a two or four channel pen recorder was used to present 
visually for monitoring all or a part of the test data as the maneuvers 
were being performed. Thus, if the control input did not produce well 
excited responses, or if a control Burface that was to be held locked 
moved inadvertently, or if the forcing control was not returned to, and 
held at, the original steady- state position this was immediately appar- 
ent. In any of these eventualities the run could be rejected immediately 
on a qualitative basis and the test repeated. 

Data Reduction 

The only flight variable requiring correction to the measured 
quantity was sideslip angle {ft ) which had to be corrected for location 
of the sideslip vane ahead of the aircraft center of gravity. The required 
correction was : 



In order to have all flight variable signals in the same angle - 
time-voltage relationship, the flight data was processed through scaling 
circuits on the analog computer prior to plotting on a pen recorder. 

The scaling factors used were determined from the calibration curves 
for each of the signal potentiometers. Where necessary, the signs of 



/? 
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the measured variables were changed to correspond to the sign of the 
corresponding output of the analog probelm. A .25 uf capacitor was 
placed across each of the output terminals of the telemeter. This 
size capacitor was found to help in eliminating high frequency noise in 
the flight data signals while at the same time having negligible effect 
on the amplitude and phase of the signals. The computer circuits 
used in correcting the sideslip angle and scaling the flight test data 
were as follows: 




In order to use the actual flight test control deflections as forcing 
functions for the analog problem it was necessary to feed them into the 
computer through a bias circuit so that the steady-state control position 
signals, which were of the order of -50 volts at the telemeter output 
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terminals, could be reduced to a zero voltage input to the analog 
problem. A potentiometer was also required in each input circuit to 
adjust the amplitude of the control deflection voltage to match the 
scaled magnitude of the control deflection on the flight test data plots. 
As in the data reduction circuits, a .25 uf capacitor was put across the 
telemeter output terminals. The above described circuits are included 
in Fig. 13. 
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DISCUSSION 

General Description of Overall Program 

In obtaining flight data, four different types of control deflections 
were used to excite the lateral modes of the NAVION. These were a 
rudder pulse, a rudder doublet, an aileron doublet, and an aileron step 
function. By means of telemeter equipment, the transient response of 
the aircraft as well as the control deflections as a function of time were 
recorded on magnetic tape at the ground station. From a large selec- 
tion of such data, one representative trace was made of each of the four 
types of forcing functions mentioned above. This was accomplished by 
processing the flight data through seeding circuits in the analog computer 
and recording the variables in time history form on a Sanborn recorder. 
A transparent overlay was then made of each of the four traces of 
scaled flight data in order to facilitate comparison of the flight and com- 
puter responses. 

Next, the lateral stability derivatives of the NAVION were 
theoretically estimated,, The lateral equations of motion based on these 
derivatives were placed on an analog computer in such a way as to enable 
independent variation of each stability derivative. The magnetically 
recorded signal corresponding to one of the four chosen control deflec- 
tions was then fed directly into the computer to force the analog. Both 
the character and magnitude of the forcing function signal were checked 
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by means of the overlay of the flight data. Once the forcing function 
was properly adjusted, the computer responses to the forcing function 
were recorded and compared to the flight responses by means of the 
same overlay. Adjustments were then made to the computer deriva- 
tives as necessary to improve the match between the flight data and 
the computer response. The procedure of forcing the computer was 
then repeated, using another one of the four forcing signals, and mak- 
ing appropriate adjustments to derivative values. This procedure 
proved to be a converging iteration process when repeated several 
times, using all four types of forcing functions. When a set of deriva- 
tives was found such that the response of the computer gave satisfactory 
matching with the flight data for all four of the forcing functions, and no 
further improvement seemed justifiable in relation to the amount of time 
involved, the value of each stability derivative was read out of the com- 
puter and tabulated. 

Prelim i nary Considerations 

Theoretical Background. At the onset, of the investigation, it 
was necessary to make a purely theoretical analysis of the proposed 
method of attack of the problem. In effect, the pertinent questions were: 
a) If a set of coefficients is determined which, when placed in the 
equations of motion, duplicate the response of the aircraft, 
can it be proved that these coefficients are unique rather than 
just one set of an entire family of possible solutions? 
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b) With two possible methods of exciting the lateral mode, i.e. , 
the rudder and the ailerons, how much data should be col- 
lected and what should be the character of the control deflec- 
tions used to excite the motion? 

Examining these questions, it was determined that one was related 
to the other. The answer to question (a) is a function of the amount of 
independent data available. Should there be a case where, in fact, an 
entire family of solutions exists, then there would be an interdependence 
between at least two of the coefficients, and one coefficient could be ex- 
pressed in terms of one or more of the other coefficients. If it were 
then possible to obtain more independent data, a point would be reached 
where there were sufficient conditions to limit the family of solutions to 
a single solution. Applying the latter case to the investigation at hand, 
it would then be known that if a set of coefficients were obtained by the 
response -matching technique, these coefficients would be unique, i.e., 
the only set of coefficients which would produce the matching of re- 
sponses. The stability derivatives conta in ed in the lumped coefficients 
would then indeed be the stability derivatives of the aircraft. The 
amount of independent data available would be in part a function of the 
control deflections used, and in addition would depend on the number of 
quantities such as roll rate, yaw rate, and sideslip angle which were 
recorded during the response. 

Fortunately, a theoretical investigation into this subject was 
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available in Ref. 4. The authors, Messrs. Klein and Sedney,, assume 
that the equations of motion are linear with constant coefficients (how- 
ever, it is pointed out that the same method of analysis can be applied 
to more general types of equations), and also that the data is free from 
experimental error. The results of their analysis is a simple working 
rule which cam be stated in words as follows: In each equation describ- 

ing a system, the number of unknown coefficients which can be deter- 
mined plus the number of known coefficients which may be checked are 
equal to the total number of measured quantities less the number of 
equations . 

This result was applied to the lateral equations of motion of the 
NAVION, with J x , J z , J xa , and known . The quantities which 
were to be recorded during the testing were roll rate, yaw rate, side- 
slip angle, and control deflections. In application of the Klein and 
Sedney criterion, roll angle, roll acceleration, yaw angle, yaw accel- 
eration, and sideslip rate are considered to be recorded in a concealed 
form in the other time histories, since they are obtainable by differen- 
tiation or integration of the appropriate trace. By following the working 
rule of the criterion, it was found that in each of the three equations of 
motion there was actually an excess of independent data under the re- 
cording conditions proposed, and therefore the solution for the coeffi- 



cients would be unique. 
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Practical Aspects. Although no theoretical limitation exists in 
the case of the lateral mode, there is a practical limitation which 
applies to the analysis of all dynamic response data. The following 
statements are extracted from Chapter 10 of Ref. 1: 

Instrument errors and unknown atmospheric disturbances 
are usually present to a sufficient degree to introduce significant 
errors into the solution of the inverse problem. As a general 
qualitative rule, the accuracy with which one can obtain each of 
the stability derivatives is associated with the importance of that 
stability derivative in influencing the motion of the aircraft dur- 
ing the test, the highest accuracy being obtained for the more 
important stability derivatives. 

Consideration of the means for improving the accuracy 
with which the less important derivatives can be obtained leads 
to these three approaches to the problem; 

1. Overall improvement of the accuracy of the recording 
instrumentation and improvement of the calibration 
and pre -flight testing techniques. 

2. Selection of test maneuvers which increase the force 
and moment contributions due to the smaller deriva- 
tives compared to the forces and moments due to the 
other derivatives. 

3. Performance of flight tests under so-called ideal 
conditions, i.e., air is smooth and stable. 

Selection of flight test instrumentation components was made 
with the requirements of (1) above in mind. As previously discussed 
under "Airborne Instrumentation,” the transducers were chosen to be 



compatible with the recording instrumentation in order to eliminate 
excessive electronic modification of the signal voltages prior to record- 
ing. This effectively reduced possible sources of error. Calibration 
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was carefully performed and the resolution and linearity of the com- 
ponents are readily apparent from the calibration curves, as presented 
in Fig. 9 and Fig. 10. Test flights were performed under conditions 
of smooth and stable air. 

The selection of suitable forcing functions for the lateral mode 
was initially considered from the viewpoint of Klein and Sedney, who 
state that theoretically, the response of the aircraft to a single forcing 
function is sufficient to establish the stability derivatives, with the 
exception of the control derivatives. Accordingly, a rudder pulse was 
selected as the forcing function most suitable for excitation of the 
lateral mode. However, the prospect of extending the method of analy- 
sis to obtain all of the control derivatives, plus the possibility of 
improving the accuracy of certain derivatives through the use of addi- 
tional forcing functions, was made very attractive when considered 
from the standpoint of paragraph (2) of the above quotation from Ref. 1. 
Four forcing functions were therefore decided upon. Considering the 
rudder pulse as the basic data, a rudder doublet was utilized to in- 
crease the amplitudes of the lateral responses, and to reduce the 
tendency of the aircraft to fall off on one wing during the transient as 
had been observed in the case of excitation by a rudder pulse. The 
aileron doublet was considered as an effective means of evaluating 

, and the aileron derivatives, whereas a time -history of 
the roll rate response to an aileron step function was considered to be 
the best means for analyzing the stability of the spiral mode. 
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.Analog of the Test Aircraft 

Theoretical Approximation of Derivatives. In order to construct 
an analog of the lateral equations of motion, it was necessary to esti- 
mate the value of the stability derivatives from whence the approximate 
value of each lumped coefficient could be determined. Anticipating the 
possibility of extremely slow convergence of the iteration process used 
on the computer, the theoretical value for each derivative was deter- 
mined as accurately as possible by several different methods and then 
a median value was taken. The results and a complete explanation of 
this work is contained in Appendix B. 

Lateral Equations of Motion. The lateral equations of motion of 
a linear aircraft are as follows: 



(c y -2 d)( 3 - a d y- + + c^Sn, = => 




For ease of construction of the analogy the equations can be 
rewritten incorporating the following changes: 

a) in terms of the highest derivative. 

b) in real time by the explicit insertion of the non-dimensional 



time constant , IT . 
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c) for simplicity, the algebraic signs changed where necessary 
so that all work from this point on could be conducted in 
terms of the absolute value of each derivative. 



The resulting equations are: 




The numerical values of the non-dimensional constants and the 
estimated derivatives were then substituted into the equations. The 
terms containing J xz were found to be so small in comparison to the 
other terms that they were immediately discarded as negligible: 

— (3 - -t- 4" - - HI 

dp = -11. & <3 +- (.14 ^ - 8-1 H -V- 1.8*1 

H' - -sn Vf- .ae><i 0 _ 3.Ub ^ -.^c, 

The constructed analog for the above equations is shown 
schematically in Fig. 11. 
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Development of Method of Analysis 

Prior to the analysis of the actual flight test data, a step by 
step procedure was developed for matching the computer response to 
flight data. First, an investigation was made to determine the effect 
of each stability derivative on the response of the system in terms of 
amplitude, period, and damping. The analog of the test airplane 
(Fig. 11) was placed on the analog computer, and the potentiometers 
set to correspond to the NAVION theoretical stability derivatives of 
Table V. An electronic function generator was used to force the 
analog with a square pulse type function, to approximate the form of a 
rudder pulse. A series of runs was then made during which the sta- 
bility derivatives were varied successively by a factor of two and then 
returned to their original theoretical values. From these tests, the 
change in the response due to the change of any one derivative was re- 
corded, and the results included in Table Ill-a. A cross-table was 
then prepared showing which derivatives should be varied to obtain a 
desired change in the response of the system. These results are 
shown in Table Ill-b. 

With the aid of the information presented in Tables Hl-a and 
in-b, an iterative procedure for the matching of computer response with 
flight data was outlined as follows: 

(a) Compare flight data and computer response curves obtained 
from the same forcing function. Note discrepancies in 



amplitude, phase, damping, etc., between corresponding 
responses. 

(b) Concentrate on the major discrepancy. By reference to 
Tables IU-a and Ill-b, determine which derivative, or 
derivatives, could be altered to correct this type of 
discrepancy. 

(c) Iteratively apply the corrections indicated in (b) above, 
commencing with the derivative having the most pronounced 
effect on the particular discrepancy. When no further 
improvement can be observed by adjusting the value of the 
first derivative, continue in a consecutive adjustment of all 
the derivatives affecting the discrepancy. 

(d) Repeat (a), (b), and (c) above for additional discrepancies 
until general agreement is obtained between flight data and 
the computer response curves. 

(e) Continue by comparing flight data and computer response 
curves obtained from a second forcing function. Again, 
repeat (a), (b), and (c) until agreement is obtained. 

(f) Repeat (e) for as many forcing functions as are available or 
practical. When the response to all forcing functions have 
been examined once, return to the first forcing function and 
repeat the entire process outlined above; continue until a 
set of derivatives is obtained which give the best possible 
response match to all forcing functions used. 
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Preliminary Test of Feasibility of Method 

A program was next performed to establish the feasibility of 
using an analog computer to curve -fit flight data and yield unique 
stability derivatives. It was considered that if the transient response 
of the computer to a rudder-pulse type of forcing function were record- 
ed, and considered as error -free "flight data',' and the computer 



potentiometers corresponding to stability derivatives then changed, 

that subsequent re -matching of the original synthetic "flight data" 

response should yield computer potentiometer settings identical to 

those of the " flight data'.' It was found that this type of analysis 

produced excellent results, with the exception of the determination of 

C , and . The responses were found to be too 

insensitive to changes in these three derivatives in order to establish 

them with satisfactory accuracy. The final value of C., contained 

71 p 

error approximately ten percent. All other derivatives were 
obtained within an error of five percent. and Cy were 

obtained with virtually no error ; this latter result was anticipated, 
from consideration of the large relative magnitude of the lumped 
coefficients containing these terms. 

Unfortunately, the above program was not extended to the use of 
an aileron input. Although the results were gratifying for the rudder 
forcing function used, they might well have been improved if the addi- 
tional forcing function data had been utilized. In particular, the use of 
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aileron forcing functions might have resulted in a better indication of 
the effect of the derivative C ^ 

The investigators considered that the results of this preliminary 
program established the feasibility of the basic method. It was also 
encouraging to note that with a minimum of practice, the iteration pro- 
cedure developed for matching the analog responses to flight data 
resulted in moderately rapid convergence. With ten stability deriva- 
tives all acting as variables in the inverse problem, a great deal of 
difficulty had been anticipated in effecting a convergence when using 
what is basically a simple iteration process, but to the degree of accur- 
acy obtained in the results mentioned above, no significant problems 
were realized. 

In future applications of this method, the procedures described in 

ii ii it 

the sections Development of Method of Analysis and Preliminary 
Test of Feasibility of Method" would be useful to an investigator from 
two standpoints : they would provide data on the sensitivity of the specific 
configuration under investigation to changes in individual derivatives, 
and would provide valuable experience in manipulating the coefficients 
to obtain the desired response. 

Application of Method to Flight Test Data 

The four control deflections which were selected and the corres- 
ponding responses of the aircraft in the three variables , (jj , and 



are shown in Figs, 15-18. As mentioned previously, a transparent 
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overlay was made of each of these four sets of data to facilitate com- 
parison of the computer response with the flight data. 

One of the important features of the instrumentation was the fact 
that signals proportional to the four control deflections were available 
for forcing the computer by playing the recording of flight data, stored 
on the magnetic tape, through the decoding component of the telemeter 
ground station. From the output of the telemeter, the circuity shown 
in Fig. 13 was used to bias the steady state signals of the controls to 
zero and to provide proper scaling of the control deflection signal 
magnitudes. Without this feature, it would have been necessary to 
obtain the forcing function for the computer from an electronic function 
generator, with limited prospects for exact duplication of the aircraft 
control deflections. The use of the exact forcing function from the flight 
data thus eliminated a possible source of error by maximizing the com- 
patibility of the data collection system and the analog computer. 

The relationship between the actual aircraft response and the 
response of the analog containing the initial theoretical derivatives, 
when each was forced by a rudder pulse, is shown in Fig. 14. With this 
initial analog response (Fig. 14) as a starting point, the attempt to 
curve -fit the analog to the flight data was commenced. 

The initial steps of the curve fitting process were the same as the 
first four parts of the procedure previously described in "Development 
of Method of Analysis Reference to Fig. 14 indicated that the most 
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obvious discrepancy was in the period of the Dutch roll oscillation. 
Since Tables IU-a and Ill-b showed C to be the derivative most 

effective in adjusting the period, the value of C was increased 

by small increments until a relatively good period match was obtained. 
The values of the rudder control derivatives (C y. , Cy, , C c ) 



were then corrected to reflect the change in C 



V* 



. The next steps 



were to attack successively the discrepancies in damping, yaw-rate 
amplitude, and roll-rate amplitude by adjusting the appropriate deriva- 
tives in accordance with the information in Tables EU-a and Ill-b. 

This procedure improved the match between the analog and the aircraft 
responses appreciably. However, in an attempt to further improve the 
matching, and at the same time gain further insight into the interplay 
of the effects of the various derivatives, all derivatives not adjusted 
in the foregoing procedures were successively altered to observe their 
effect on the analog responses. When no change, or an adverse change, 
in the responses was encountered, the derivative under investigation 
was returned to its original value. Those changes which appeared 
beneficial were retained. In each case, the effect of altering a deriva- 
tive, and the extent of the change were recorded. 

After a reasonably close match between the analog and aircraft 



responses to a rudder pulse was obtained, the analog was next forced 
by a rudder doublet, with the computer potentiometers still set on the 
derivatives established by matching responses to the rudder pulse. 
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The procedure which was then followed to minimize the discrepancies 
between the analog and aircraft responses was identical to that used 
for response to a rudder pulse, and again produced an improvement in 
the analog responses. 

The responses to an aileron doublet and aileron step function 
were then investigated in a like manner; that is, by first concentrating 
on the derivative changes indicated in Tables ffi-a and Ill-b, and then 
investigating the individual effects of all other derivatives. During the 
initial analysis of the aileron data, it was determined that an error 
existed in the aileron signal due to cable stretch in the control system. 
Corrections required for this error are discussed in Appendix C. 

During the conduct of the actual flight tests, it had been noted 
that in response to an aileron step function, the aircraft developed a 
nose-down attitude within two seconds after the application of the 
aileron displacement. This was caused by the large roll angle that 
developed. The resulting errors in the /& and traces were 

readily apparent and consequently, these two variables were disregarded 
in this case. The only matching that was attempted with aileron step 
forcing functions were the amplitude of the initial peak of the roll rate 
response, and the neutral stability of the spiral mode as indicated by 
the constant mean value of roll rate after the roll rate response had 
reached its initial peak value. The rapid development of the nose-down 
attitude when the aircraft was forced by an aileron step function could 
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have been avoided by the use of an initial condition, such as a steady 
sideslip. This would have extended the time during which the airplane 
was in the linear approximation range. Since the equations of motion 
of the analog are developed for perturbations from steady state, a 
similar initial condition would not be required when forcing the 
computer. 

After completion of one series of iterations for the values of the 
stability derivatives, the procedure was repeated for all forcing func- 
tions, again starting with the rudder pulse. It was noted that good 
matching was quickly obtained for any one of the forcing functions 
taken individually, but a shift to the other forcing functions offered a 
relatively poor match in the early stages. Gradually, however, as a 
result of the experience obtained, and by reference to the records of 
the effects of varying the individual derivatives, new procedures for 
establishing the value of some of the derivatives became apparent. 
These procedures represented a more definite method of evaluating 
the corresponding derivatives than did the influences indicated in 
Tables ni-a and Ill-b, because in each case a particular response 
characteristic was controlled by an individual derivative. The deriva- 
tives and their individual effects are as follows: 

C n ^ - The period of the Dutch roll mode is controlled by 
r' 

this derivative. This derivative has more than ten 
times the effect on the period than does any of the 



other derivatives. 
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C n has a controlling effect on the damping of the 
r 

Dutch roll mode when the analog is forced by a 
rudder deflection. 

Changing this derivative changes the effectiveness of 

the rudder in producing yaw. Therefore, Cn * 

®r 

can be established by matching /S and (j/ ampli- 
tudes due to a rudder deflection. 

The magnitude of this derivative controls the match- 
ing of response in two respects; phase and ampli- 

tude. Phase error in the initial (fi response of the 
analog is affected considerably by this derivative. 

When the airplane is forced with an aileron step func- 
tion, the Dutch roll mode is excited to such a small 
degree that the initial peak of the <p response can 
be attributed to the rolling mode and the magnitude 
of Cw . For a given Ci , . can be 

determined by matching the initial peak of the (p 
response. Determination of the precise value of 
Cl is explained below. 

Just as C_ controls the damping of the Dutch roll 

T - 

mode, C^ controls the damping of the rolling mode. 

With C\ , adjusted as 

°<x. 

be determined by matching the response to an 



described above, C]^ can 
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aileron doublet, especially during the latter half 
of the aileron input. 

The effect of this derivative is readily apparent in 

the 0 response as an initial adverse 0 due 

to a rudder deflection. Matching the initial form of 

the (ft response will thus establish Ci^» 

In the response of the computer analog to a rudder 

deflection, Cnp has only a small effect. However, 

an upper and lower limit can be established for Cn^ » 

beyond which the character of the ^ trace is 

noticeably altered in magnitude, in the case of an 

aileron deflection, the effects of C n and C n * 

P 

are so similar that one can be replaced by the other 
to a large extent. If an iteration is performed between 
these two derivatives, however, maintaining C n ^ 
within the limits developed during a rudder deflection, 
until the best match is obtained in due to an 

aileron doublet, and C n ^ can be determined 

within reasonable limits of accuracy. 

Only a very rough value of this derivative can be 
established. Gross changes in will effect the 

stability of the spiral mode. The 0 response to 
an aileron step function is an indication of the spiral 
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mode stability: a constant mean value of (p 
indicates neutral stability, a decreasing mean value 
indicates positive stability and an increasing mean 
value indicates negative stability of the spiral mode. 
Since the other three derivatives which affect the 
spiral mode ( Cn^ , Cn^ . , C]^ ) can be determined 

by the methods outlined above, the value of can 

be adjusted to match the character of the mean 
response to an aileron step function. 

Because of the insensitivity of the responses to 
changes in these derivatives, no refinement of the 
theoretically determined values was possible. 



The above relations are summarized in Table IV. 



By the time sufficient information had been accumulated to 
develop the procedures outlined above, the initial discrepancies between 
the analog and the aircraft responses, for all four forcing functions, 
had been appreciably reduced. Subsequent use of the iteration proced- 
ure, in the case of each forcing function and its related responses, was 
confined primarily to adjustments in the derivatives which were based 
on the relationships of Table IV. In addition, when working with the 
responses to the rudder forcing functions, care was taken not to make 
large changes in those derivatives which had been established primarily 
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by matching responses to aileron forcing functions. This restriction 
was found necessary, since alteration of the aforementioned derivatives 
generally resulted in only minor improvements in the rudder cases, but 
had an appreciable adverse affect on the matching of the responses in 
the aileron cases. The same restriction was applied when matching 
responses to aileron forcing functions ; in this case, no extensive 
changes were made in those derivatives which had been established 
primarily by matching responses to rudder forcing functions. 

When a set of derivatives was found such that the responses of 
the computer gave satisfactory matching with the flight data for all 
four of the forcing functions, and no further improvement seemed justi- 
fiable in relation to the amount of time it would involve, the final value 
of each derivative was read off of the computer potentiometers and 



recorded. 
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RESULTS 



The final values of the lateral stability derivatives of the NAVION 
as determined by the method of this report are tabulated and compared 
with the initial theoretical values in Table V. The matching of the 
responses of the final computer configuration to the flight test data for 
all four types of forcing functions are presented in Figs. 15-18. 

During the iteration process, a better match had been obtained to 
the responses for each of the four forcing functions taken individually, 
but no single set of derivatives gave as good a match to the responses 
of all four forcing functions taken as a group as did the final set of 
derivatives listed in Table V. The major discrepancies remaining in 
the final matches are the generally low response of the analog in the 
/* variable for all forcing functions ; the too large amplitude re- 
sponse in the ijj variable for rudder doublet forcing functions ; and 
the initial (j/ response to the aileron doublet. It is felt that had 
more time been available, further investigation into the interrelated 

effects of the cross derivatives C„ and Ci would have resulted 

Up A r 

• 9 

in nearly perfect response matches in the (J/ and (ft variables 
when using the rudder forcing functions. 

The accuracy of the final value of each derivative is directly 
related to the sensitivity of the system response to changes in the 
derivative. The derivatives Ci » Ci , , C_ , C n . , Ci , 

U V h * 
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Cl and Cn are judged to be very accurately determined by the 
ip > n r 

method. (See Appendix C concerning the accuracy of the specific 
values for Cl, and Cq of this report.) The values for Cl » 

« <)<x, 1 

Cn and C are susceptible to somewhat larger error due to 

V w 

the smaller effects they produce on the response, and also due to 
overlapping effects in the case of C n and Cn t • Finally, the 

P 0 CL, 

value of the derivatives Cy and Cy could not be established 

V 0 A/ 

with satisfactory accuracy. Other types of flight maneuvers might be 
advantageous in establishing these latter derivatives. 

The proposed method of extracting stability derivatives from 
transient flight data was determined to be completely feasible. As 
mentioned previously, there are inherent difficulties in the process of 
extracting all of the derivatives from flight data regardless of the 
method used. For certain applications, however, no better criterion 
could be established for obtaining the values of the stability deriva- 
tives than to match the response of the analog containing these deriva- 
tives to the response of the actual aircraft. In the field of servomech- 
anisms, for example, this criterion would be considered to be of 
prime importance. 

Another advantage of the method would be its economy in time, 
both in obtaining flight data and in extracting the derivatives from that 
data. In obtaining the flight data, any arbitrary forcing function can be 
used, providing it offers sufficient excitation of the modes of aircraft 
motion. Use of the telemeter data collection system permitted 
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convenient monitoring of flight data during the flights, and thus when 
necessary, instructions could be given to the pilot, to improve the 
responses obtained. Normally this would allow the recording of all 
the data for one flight condition in a matter of minutes. The ease of 
extracting the derivatives would tend to increase with practice. Un- 
doubtedly, Table IV, which contains procedures for establishing 
several of the derivatives, could be expanded with the aid of more 
time and experience with the method. For this reason, it is felt that 
better matching than is illustrated in this report would definitely be 
possible. It is important to note that the results of this investigation 
were obtained for the aircraft in the cruise condition only. Because of 
the large variation of certain stability derivatives with changes in 
flight conditions, further investigation of the method throughout a 
range of flight conditions is recommended. There is virtually no 
chance of computational errors affecting the results, for once the analog 
has been completed the computations are performed by the computer. 

Although the theoretical estimates of the derivatives were made 
as accurately as possible in this investigation, it was found that the 
rate of convergence of the iteration process was sufficiently rapid in 
the first stages to suggest that very rough estimates of the derivative 
values would be sufficient for use as initial computer settings. The 
time involved in refining the rough estimates would more than likely 
exceed the amount of computer time necessary to accomplish the same 
degree of refinement. 
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Certain limitations of the method axe readily apparent. The 
exact form of the equations of motion must be known in order to con- 
struct the analog. The application of the method to more complex 
aircraft would require further investigation and refinement; for ex- 
ample, any nonlinear effects of the aircraft system would have to be 
duplicated by electronic circuitry in the analog. The accuracy obtain- 
able for those derivatives to which the response is insensitive is 
undoubtedly of low order, but the possibility does exist of using differ- 
ent types of flight tests to establish these values. If these tests 
were compatible with the proposed data collection and reduction 
method, the additional effort could be minimized. 
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CONCLUSIONS AND RECOMMENDATIONS 



The proposed method of extracting stability derivatives from 
transient flight data was proved to be feasible within the scope of this 
investigation. Matching of responses between the test aircraft and the 
equations of motion containing the derivatives is an excellent criterion 
for many applications of the derivatives. Due to the compatibility of 
the data collection system and the analog computer, the method is rela- 
tively fast and free from computational error. Only rough approxima- 
tions to the values of the stability derivatives are needed to construct 
the analog. In general the accuracy with which each of the derivatives 
is obtained is associated with the importance of that stability derivative 
in influencing the motion of the aircraft during the test. The use of 
four different types of control deflections as forcing functions is con- 
sidered a minimum requirement for accurate determination of all 
lateral derivatives. Subject to these restrictions, the exact form of 
each particular control deflection is arbitrary. Application of this 
method should include the preliminary procedures described under 
"Development of Method of Analysis " and " Preliminary Test of 
Feasibility of Method. " 

Application of this method to non-linear aircraft would require 
further investigation and refinement. Any non-linear effects in the 
aircraft system would have to be duplicated by electronic circuitry in 
the analog. 
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As a result of the investigation, the following recommendations 
are proposed: 

1. Investigate the method at other flight conditions and with more 
complex aircraft. 

2. Investigate the use of other control deflections to obtain 
responses which are more strongly influenced by the less 
sensitive derivatives. As an example, Cl^ might be more 
accurately obtained from an analysis of the response to a 
rudder step function. 

3. Apply an initial condition such as a steady sideslip, prior to 
forcing the airplane with an aileron step function. 

4. Augment the data (Table III) which indicates the effect(s) of 
the individual derivatives on the response of the system by 
conducting a similar variation of the derivatives when forcing 
the computer problem with an aileron function. 

5. Improve the test aircraft aileron instrumentation by recording 
both left and right aileron angles and converting to average 
aileron angle by means of analog computer scaling circuits. 
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TABUS II 



PHYSICAL CHARACTERISTICS CF NAVION 



A. WING 



1. 


Total Wing Area, Sy (includes flaps, ailerons 
and 19o87 Ft3* covered by the fuselage) 


1 8b .337 Ft. 2 


2. 


Span by 


33o378 Ft. 


3. 


MAC 


68.352 in. 


1*. 


Angle of Incidence 
Root, i r 


♦2° 




Tip, i t 


-1° 


5. 


Twist 

Aerodynamic 

Geometric 


2° 31' 
3° 


6. 


Airfoil Section 
Root 
Tip 


NACA Utl5R 
NACA 61*1CR 


7<> 


Aspect Ratio, ARy 


6.014* 


8o 


Taper Ratio, A w 


0.526? 


9. 


Dihedral 


7o5° 


10. 


Root Chord 


7.2 Ft. 


iio 


Tip Chord 


3.92 Ft. 


AU£RCN (one Aileron only) 


120°-120° 


(Beaded skin, Frlse Nose Balance, 
No Trailing Edge Extension) 


Wheel Throw 
Fixed bend tab 
right aileron. 


1. 


Area, S a 


2.161 Ft. 2 


2o 


Span, b a 


61.987 in. 


3. 


Ailarcn Deflection, ^ a 


30° DP, 20° DN 


ho 


Boost 


None 



5. 



Trla Tab 



Fixed bend tab on 
right aileron. 



